Directed Flow of Baryons in Heavy-Ion Collisions 
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The collective motion of nucleons from high-energy heavy-ion collisions is analyzed within a rel- 
ativistic two-fluid model for different equations of state (EoS). As function of beam energy the 
theoretical slope parameter F y of the differential directed flow is in good agreement with experi- 
mental data, when calculated for the QCD-consistent EoS described by the statistical mixed-phase 
model. Within this model, which takes the deconfinement phase transition into account, the exci- 
tation function of the directed flow (P x ) turns out to be a smooth function in the whole range from 
SIS till SPS energies. This function is close to that for pure hadronic EoS and exhibits no minimum 
predicted earlier for a two-phase bag-model EoS. Attention is also called to a possible formation of 
nucleon antiflow (F y < 0) at energies > 100 A-GeV. 

PACS numbers: 24.85. +p, 12.38.Aw, 12.38Mh, 21.65. +f, 64.60.-i, 27.75. +r 



I. INTRODUCTION 

Collective flows of various types (radial, directed, el- 
liptic,...) observed experimentally in heavy-ion collisions 
reveal a space-momentum correlated motion of strongly 
interacting nuclear matter. This collective motion is es- 
sentially caused by the pressure gradients arising during 
the time evolution in the collision, and hence opens a 
promising way for obtaining information on the equation 
of state (EoS) and, in particular, on a possible phase 
transition. Recently, this feature has stimulated a large 
number of experimental and theoretical investigations on 
flow effects (cf. review articles Jl|,||). 

Manifestations of the deconfinement phase transition 
have been considered already some time ago by Shuryak, 
Zhirov D and van Hove M. Since a phase transition 
slows down the time evolution of the system due to soft- 
ening of the EoS, the authors expect around some crit- 
ical incident energy a remarkable loss of correlation be- 
tween the observed particle momenta and the reaction 
plane, and hence a reduction of the directed flow. As- 
suming a first-order phase transition Hung, Shuryak ||] 
and Rischke et al. || have recently obtained quantitative 
predictions for heavy-ion collisions. For an expanding 
fireball Hung and Shuryak expect the softest point ef- 
fect around Ei a b — 30 A-GeV. Within a one-fluid hydro- 
dynamic model Rischke et al. show that the excitation 
function of the directed flow exhibits a deep minimum 
near Ei a b = 6 A-GeV. However, preliminary experimen- 
tal results in this energy range do not confirm these 
predictions. In the following, we report on a study of the 
directed flow within a two-fluid hydrodynamic model ||] 
for the statistical mixed-phase EoS ^,[l(| which is ad- 
justed to available lattice QCD data. 



II. EQUATION OF STATE WITHIN THE 
MIXED-PHASE MODEL 

Our consideration is essentially based on the recently 
proposed Mixed-Phase (MP) model |p|, |l0f which is consis- 
tent with available QCD lattice data [ |ll| . The underly- 
ing assumption of the MP model is that unbound quarks 
and gluons may coexist with hadrons forming a homoge- 
neous quark/gluon-hadron phase. Since the mean dis- 
tance between hadrons and quarks/gluons in this mixed 
phase may be of the same order as that between hadrons, 
the interaction between all these constituents (unbound 
quarks/gluons and hadrons) plays an important role and 
defines the order of the phase transition. 

Within the MP model [g,|l0| the effective Hamilto- 
nian is expressed in the quasiparticle approximation with 
density-dependent mean-field interactions. Under quite 
general requirements of confinement for color charges, 
the mean-field potential of quarks and gluons is approx- 
imated by 



U q (p) = U g {p) = ^ ; 7 > 
with the total density of quarks and gluons 



(1) 



P = Pq + Pg + 



E 



iPj 



where p q and p g are the densities of unbound quarks 
and gluons outside of hadrons, while pj is the density of 
hadron type j and Uj is the number of valence quarks 
inside. The presence of the total density p in (Q) implies 
interactions between all components of the mixed phase. 
The approximation (Q) mirrors two important limits of 
the QCD interaction. For p — ► 0, the interaction poten- 
tial approaches infinity, i.e. an infinite energy is neces- 
sary to create an isolated quark or gluon, which simulates 
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the confinement of color objects. In the other extreme 
case of large energy density corresponding to p — > oo, we 
have U q = U g = which is consistent with asymptotic 
freedom. 

The use of the density-dependent potential ([[]) for 
quarks and the hadronic potential, described by a mod- 
ified non-linear mean- field model (l|, requires certain 
constraints to be fulfilled, which are related to ther- 
modynamic consistency |^,[lCj. For the chosen form of 
the Hamiltonian these conditions require that U g (p) and 
U q (p) do not depend on temperature. From these condi- 
tions one also obtains a form for the quark-hadron po- 
tential ||. 

A detailed study of the pure gluonic SU (3) case with 
a first-order phase transition allows one to fix the values 
of the parameters as 7 = 0.62 and A 1/(37+1) = 250 MeV. 
These values are then used for the SU (3) system includ- 
ing quarks. As is shown in Fig.l for the case of quarks 
of two light flavors at zero baryon density (ng = 0), the 
MP model is consistent with lattice QCD data providing 
a continuous phase transition of the cross-over type with 
a deconfinement temperature T^ ec = 153 MeV. For a 
two-phase approach based on the bag model a first-order 
deconfinement phase transition occurs with a sharp jump 
in energy density e at Td ec close to the value obtained 
from lattice QCD. 
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FIG. 1. The reduced energy density e/ess and pressure 
p/psb (the ess and psb are corresponding Stephan-Boltz- 
mann quantities) of the SU(3) system with two light flavors 
for n,B = calculated within the MP (solid lines) and bag 
(dashed lines) models. Circles and squares are lattice QCD 
data obtained within the Wilson |L3| and Kogut-Susskind |ll|] 
schemes, respectively. 

Though at a glimpse the temperature dependences of 
the energy density e and pressure p for the different ap- 
proaches presented in Fig.l look quite similar, there are 
large differences revealed when p/e is plotted versus e (cf. 
Fig. 2, left panel). The lattice QCD data differ at low e, 
which is due to difficulties within the Kogut-Susskind 
scheme [[l4| in treating the hadronic sector. A particular 
feature in the MP model is that, for ns = 0, the softest 
point of the EoS, defined as a minimum of the function 



p( £ )/ £ @j is n °t ver y pronounced and located at com- 
paratively low values of the energy density: esp ~ 0.45 
GeV/fm 3 , which roughly agrees with the lattice QCD 
value [|l3| . This value of e is close to the energy density 
inside the nucleon, and hence, reaching this value indi- 
cates that we are dealing with a single big hadron con- 
sisting of deconfined matter. In contradistinction, the 
bag-model EoS exhibits a very pronounced softest point 
at large energy density esp ~ 1.5 GeV/fm 3 ||||. 
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FIG. 2. The (e,p/e)-representation of the EoS for the 
two-flavor 5(7(3) system at various baryon densities hb- No- 
tation of data points and lines is the same as in Fig.l. 

The MP model can be extended to baryon-rich sys- 
tems in a parameter- free way |P,|Io|| . As demonstrated in 
Fig. 2 (right panel), the softest point for baryonic matter 
is gradually washed out with increasing baryon density 
and vanishes for ub > 0.4 rig (no is normal nuclear mat- 
ter density) . This behavior differs drastically from that of 
the two-phase bag-model EoS, where esp is only weakly 
dependent on ng [§|,D- It is of interest to note that the 
interacting hadron gas model has no softest point at all 
and, in this respect, its thermodynamic behavior is close 
to that of the MP model at high energy densities [ p"o| . 

These differences between the various models should 
manifest themselves in the dynamics discussed below. 



III. TWO-FLUID HYDRODYNAMIC MODEL 

In contrast to the one-fluid hydrodynamic model, 
where local instantaneous stopping of projectile and tar- 
get matter is assumed, a specific feature of the dynami- 
cal two-fluid description is a finite stopping power. Ex- 
perimental rapidity distributions in nucleus-nucleus colli- 
sions support this specific feature of the two-fluid model. 
In accordance with the total baryonic current and 
energy-momentum tensor are written as 



(2) 
(3) 
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where the baryonic current J£ — n a u^ and energy- 
momentum tensor T£" of the fluid a are initially asso- 
ciated with either target (a = t) or projectile (a = p) 
nucleons. Later on - while heated up - these fluids con- 
tain all hadronic and quark-gluon species, depending on 
the model used for describing the fluids. The twelve inde- 
pendent quantities (the baryon densities n a , 4- velocities 
u£ normalized as u afl u^ = 1, as well as temperatures 
and pressures of the fluids) are obtained by solving the 
following set of equations of two-fluid hydrodynamics Iq] 



« = 

U ^ ± a — ot 



where the coupling term 



n p n t (Vrel J da NN ^ NX {s) {jp-p a ) u \ 



(4) 
(5) 



(0) 



characterizes friction between the counter-streaming flu- 
ids. Here, n 8 a and {p — p a ) denote respectively the scalar 
density of the fluid and the 4-momentum transfer gained 
by a particle of the fluid a after collision with a par- 
ticle of the counter-streaming fluid. The cross sections 
d&NN—>NX take into account all elastic and inelastic in- 
teractions between the constituents of different fluids at 
the invariant collision energy s 1 / 2 with the local relative 
velocity V re i — [s(s — Am 2 N )] 1 / 2 /2m 2 N . The average in (||) 
is taken over all particles in the two fluids which are as- 
sumed to be in local equilibrium intrinsically . The set 
of Eqs. (|]) and (||) is closed by an EoS, which is naturally 
the same for both colliding fluids. 

The friction term F£ in Eq. (|J) originates from both 
elastic and inelastic NN collisions. The latter give rise 
to a direct emission of mesons in addition to the thermal 
mesons in the fluids. In the present version only for the 
pions the direct emission is included by the additional 
equations 



dvK = n s p n s t 



V r , 



V r , 



drr 



NN- 



tX 



da 



NN^-rrX P„ 



(7) 



where p„ is the 4-momentum of an emitted direct pion. 
These equations together with (||) provide the total 
energy-momentum conservation 



V 



TD = o 



(9) 



It is assumed []S| that in the subsequent evolution these 
direct pions interact neither with the fluids nor with each 
other. This is a reasonable assumption at relativistic 
energies, simulating a long formation time of these direct 
pions. At moderate energies, where the latter argument 
does not hold in general, the number of direct pions is 
negligible compared to the number of thermal pions. 

For the calculation of the friction force (||) , approxima- 
tions of N-N cross-sections are used. It was found fllal 



that a part of the friction term, which is related to the 
transport cross-section, may be well parametrized by an 
effective deceleration length A e g with a constant value 
A ff ~ 5 fm. However, there are reasons to consider A e ff 
as a phenomenological parameter, as it was pointed out 
in ||. First, the value of A e gf is highly sensitive to the 
precise form of parameterization of the free cross-sections 
which, in addition, may be essentially modified by in- 
medium effects. Furthermore, the model neglects the in- 
teractions of direct pions both with each other and with 
baryons, as well as the direct emission of other mesons 
which are produced quite abundantly at SPS energies. 
Due to all these effects the stopping power at SPS ener- 
gies is somewhat underestimated [p|. This shortcoming 
of the model is cured by an appropriate choice of the A c g 
value as 

A c ff = a exp(—by/s) 

with a = 6.6 fm and b = 0.106 GeV -1 adjusted to the 
rapidity distributions of nucleons and pions in central 
Au+Au collisions at AGS and SPS energies. 

Following the original paper g), it is assumed that a 
fluid element decouples from the hydrodynamic regime, 
when its baryon density ub and densities in the eight 
surrounding cells become smaller than a fixed value nj. 
A value n/ = 0.8no was used for this local freeze-out 
density which corresponds to the actual density of the 
frozen-out fluid element of about 0.5no to 0.7uq. 



IV. COLLECTIVE FLOW FROM HEAVY-ION 
COLLISIONS 

For central nucleus-nucleus collisions only the isotropic 
transverse expansion, or transverse radial flow, develops 
due to the azimuthal symmetry of a system. The pres- 
ence of the reaction plane for non-central collisions de- 
stroys this symmetry and gives rise to various patterns 
of collective motion generated by compressed and excited 
nuclear matter created during the collision. For example, 
the directed (or sideward) flow characterizes the deflec- 
tion of emitted hadrons away from the beam axis within 
the reaction plane. In particular, one defines the dif- 
ferential directed flow by the mean in-plane component 
(p x (y)) of the transverse momentum at a given rapidity 
y. This deflection is believed to be quite sensitive to the 
elasticity or softness of the EoS. 

The (px(y)) distributions of baryons are shown in Fig. 3 
for Au+Au collisions at Ei a b — 10 A-GeV calculated for 
different EoS at an impact parameter 6 = 3 fm. In gen- 
eral, the characteristic S'-shape of the distribution is re- 
produced, demonstrating a definite anti-correlation be- 
tween nucleons bounced-off from the target and projec- 
tile regions. One should keep in mind that the protons 
bound in observed complex particles (e.g., in deuterons) 
are not excluded in our calculations. Therefore, all hy- 
drodynamic results should be compared to the triangle 
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points in Fig. 3, where nucleons from complex particles 
do contribute. The MP and interacting hadronic mod- 
elsQ give similar results, both getting into error bars of 
these triangle points, though the flow in the MP model is 
slightly lower due to softening of EoS near the crossover 
phase transition. This softening is stronger for the bag- 
model EoS. However, one should note that different ver- 
sions of 
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FIG. 3. Differential directed flow of nucleons in the re- 
action plane as a function of rapidity in semi-central (the 
trigger transverse energy Et ~ (200 — 230) GeV) Au + Au 
collisions at the energy 10 A-GeV. Three curves are calculated 
within relativistic two-fluid hydrodynamics for an impact pa- 
rameter b = 3 fm and different EoS: for the MP model (solid 
line), interacting hadron gas model (dashed) and two-phase 
bag model (dot-dashed). Circles are experimental points for 
identified protons, triangles correspond to a nucleon flow es- 
timate based on the measurement of Et and the number of 
charged particles No Experimental points marked by 

full symbols are measured directly, open ones are obtained by 
reflecting at the mid-rapidity point. 

transport codes, which do not account for a phase tran- 
sition, also give a reasonable description of (p x (y)) i e -9-, 
see the comparison with RQMD results in [[16]]). There- 
fore, convincing evidence on a possible phase transition, 



*The interaction in the hadron model is taken into account 
in the same manner as that in the hadronic sector of the MP 
model iRIlC 



based solely on the data at a single bombarding energy, 
is hardly possible. 

The rapidity dependence of the mean in-plane trans- 
verse momentum can be quantified in terms of the deriva- 
tive at mid-rapidity 



Fy = 



d (px(y)) 



dy 



(10) 



which is quite suitable for analyzing the flow excitation 
function. The slope parameter F y calculated with the 
MP-model EoS is presented in Fig. 4 (upper panel) to- 
gether with available experimental points covering the 
whole range of incident energies. The results describe 
correctly the decrease of F y with increasing energy and 
shows essentially no dependence on impact parameter 
for semi-central collisions. It was shown experimen- 



tally 1 17 1 that the directed flow is larger for heavier frag- 
ments. As mentioned before, hydrodynamic calculations 
deal with primordial nucleons, and hence they describe 
the mean value of F y for free nucleons and nucleons 
bound in deuterons and heavier fragments. Therefore, 
our hydrodynamic results lie between the experimental 
points for identified protons (open circles in Fig. 4) and 
the data [17) (full circles in Fig. 4) for intermediate mass 
fragments 17 This effect is particularly strong for ener- 
gies below Eiab ~ 1 A-GeV, where the baryonic flow is 
largest (and heavy fragments in the mid-rapidity range 
are abundant). 

To our knowledge no other hydrodynamic calculations 
of Fy(Ei a b) have been reported. Therefore, we compare 
our results with transport calculations in the lower panel 
of Fig. 4. The ARC and ART are cascade models, while 
the RQMD takes also into account mean-field effects. 
Though all these models agree with experimental data 
at Eiai, sa 10 A-GeV (considered as a reference point), 
values of F y at lower energies are clearly underestimated, 
as is evident from comparison with preliminary results of 
the E895 Collaboration (see empty squares in Fig. 4). 
Recently, a good description of experimental points (in- 
cluding the E895 data) was reported within a relativistic 
BUU (RBUU) model [^l| . The good agreement with ex- 
periment was achieved by a special fine tuning of the 
mean fields involved in the particle propagation. 



^For complex particles the value F y was deduced from p x per 
baryon. The FOPI data for intermediate mass fragments are 
often scaled by factor 0.7 to make absolute values comparable 
to those for p, d, a jl7j . We do not use this scaling factor. Note 
that everywhere we deal with F y defined by Eq.(|lo|), which 
differs from the frequently used flow parameter F = F y y cm - 
The extra beam-rapidity factor y cm obscures a relative role of 
dynamical effects in the energy dependence of the directed- 
flow slope parameter. 
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FIG. 4. Excitation function of the slope parameter F y for 
baryons from Au + Au collisions within two-fluid hydrody- 
namics for the MP EoS (upper panel) and within different 
transport simulations (lower panel). Open symbols are ex- 
perimental points for identified protons (see data collection 
in P Jl^Jl^ l), filled circles, triangles and squares correspond 
to the flow parameter measured for intermediate mass frag- 
ments jl7| and for light particles p, d, a [ p"9pc| ] . The results 
of transport calculations for three different codes are given 
by the thin solid (RQMD), dashed (ARC) and dot-dashed 
(ART) lines (cited according to @) . The solid line (RBUU) 
is taken from 

It is of interest to mention that the calculated value 
of F y for the baryon flow becomes negative (antiflow) for 
beam energies > 100 A-GeV, while the experiment p^ ] 
gives a small but positive value even at 158 A-GeV. The 
reason of this antiflow is a wiggle in the (p x (y)} distri- 
bution arising in hydrodynamic results within a narrow 
mid-rapidity interval \6y\ < 1 due to a peculiar interplay 
between the transverse radial and directed flows. The 
possibility of such an effect was noticed in [ p3[ some time 
ago and later also observed in the UrQMD transport cal- 
culations p4j . However, actual measurements have been 
taken at larger rapidities and then extrapolated into this 
unmeasured region |22]. Therefore, more accurate data 
in the mid-rapidity region are necessary to clarify this 



behavior. 

The directed flow can be characterized by another 
quantity which is less sensitive to possible rapidity fluc- 
tuations of the in-plane momentum. Such a quantity is 
the average directed flow which is defined by 



(P*) 



dp x dp y dy p x E 



dp x dp y dy E 



d 3 N 



'I I' 3 ) 



(11) 



where the integration in the cm. system runs over the ra- 
pidity region [0, y C m}- The calculated excitation functions 
for the average directed flow of baryons within different 
models are shown in Fig. 5. Conventional (one- fluid) hy- 
drodynamics for pure hadronic matter B results 
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FIG. 5. The excitation function of the average directed flow 
for baryons from Au + Au collisions. Two-fluid hydrody- 
namics with the MP EoS at the impact parameter 3 fm is 
compared with the corresponding results of one-fluid |j| (up- 
per panel) and three-fluid (lower panel) hydrodynamics 
with the bag-model EoS. One-fluid calculations both with and 
without the phase transition (PT) are displayed. 

in a very large directed flow due to the inherent instan- 
taneous stopping of the colliding matter. This instanta- 
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neous stopping is unrealistic at high beam energies. If 
the deconfinement phase transition, based on the bag- 
model EoS |@] , is included into this model, the excitation 
function of (P x ) exhibits a deep minimum near Ei a f, w 6 
A-GeV, which manifests the softest-point effect of the 
bag-model EoS depictured in the right panel of Fig. 2. 

The result of two-fluid hydrodynamics with the MP 
EoS noticeably differs from the one-fluid calculations. Af- 
ter a maximum around 1 A-GeV, the average directed 
flow decreases slowly and smoothly. This difference is 
caused by two reasons. First, as follows from Fig. 2, the 
softest point of the MP EoS is washed out for ns > 0.4. 
The second reason is dynamical: the finite stopping 
power and direct pion emission change the evolution pat- 
tern. The latter point is confirmed by comparison to 
three-fluid calculations with the bag EoS |^5| plotted in 
the lower panel of Fig. 5. The third pionic fluid in this 
model is assumed to interact only with itself neglecting 
the interaction with baryonic fluids. Therefore, with re- 
gard to the baryonic component, this three-fluid hydro- 
dynamics |^5|,^6| is completely equivalent to our two-fluid 
model and the main difference is due to the different EoS. 
As seen in Fig. 5, the minimum of the directed flow excita- 
tion function, predicted by the one-fluid hydrodynamics 
with the bag EoS, survives in the three-fluid (nonunified) 
regime but its value decreases and its position shifts to 
higher energies. If one applies the unification procedure 
of [^H, which favors fusion of two fluids into a single 
one, and thus making stopping larger, three-fluid hydro- 
dynamics practically reproduces the one-fluid result and 
predicts in addition a bump at Ei a i, « 40 A-GeV. 



V. CONCLUSIONS 

We have studied relativistic nuclear collisions within 
3D two-fluid hydrodynamics combined with different 
EoS, including that of the statistical mixed-phase 
model of the deconfinement phase transition, developed 
in |],[T(|. It has been shown that the directed flow ex- 
citation functions F y and (P x ) for baryons are sensitive 
to the EoS, but this sensitivity is significantly masked by 
nonequilibrium dynamics of nuclear collisions. Neverthe- 
less, the results indicate that the widely used two-phase 
EoS, based on the bag model |||| and giving rise to a 
first-order phase transition, seems to be inappropriate. 
The neglect of interactions near the deconfinement tem- 
perature results in an unrealistically strong softest-point 
effect within this two-phase EoS. In fact, its prediction of 
a minimum in (P x ) {Ei a b) near E\ a \, w 6 A-GeV has not 
been confirmed experimentally [Q. However, accurate ex- 
perimental investigations of the differential directed flow 
and flow excitation functions in the energy region be- 
tween AGS and SPS are still highly demanded not only 
in searching for a shifted minimum of (P x ) (Ei a b), but 
also in clarifying the physics of a possible negative slope 
(antiflow) of the baryonic directed flow F y . This antiflow 



is particularly sensitive to the EoS. While for the EoS in 
the MP model the antiflow is predicted at incident ener- 
gies only above 100 A-GeV, it occurs already at 8 A-GeV, 
when the bag EoS is used |p5|] . 

In this respect a dramatic phenomenon of the cracked 
nut proposed recently as a hydrodynamic signature of the 
QCD phase transition at RHIC and LHC energies |27 28 
looks questionable. The authors argue that the softest 
point in the EoS may lead to the development of two 
shells at the edge of the almond- like initial fireball, which 
are then cracked by internal pressure and separated, re- 
sulting in a specific flow pattern. However, this specula- 
tion was based on the bag-model EoS. The application of 
the QCD-consistent EoS of the MP model to this prob- 
lem would be interesting. 

The directed flow is the first coefficient in the Fourier 
decomposition of the azimuthal momentum distribution 
of particles ^9|. The second coefficient, the elliptic flow, 
is expected to be more sensitive to the EoS and some 
hints of the phase transition have been indicated by an 
analysis of the measured excitation function for the ellip- 
tic flow (see review articles KL0|). The study of the ellip- 
tic flow within two-fluid hydrodynamics for the mixed- 
phase model EoS is in progress. 



ACKNOWLEDGMENTS 

We are grateful to V. Russkikh for making his hydro- 
dynamic code available to us. We thank P. Danielewicz, 
B. Friman and E. Kolomeitsev for useful discussions. 
Yu.B.L, E.G.N, and V.D.T. gratefully acknowledge the 
hospitality at the Theory Group of GSI, where this work 
has been done. This work was supported in part by DFG 
(project 436 RUS 113/558/0) and RFBR (grant 00-02- 
04012). Yu.B.L was partially supported by RFBR grant 
00-15-96590. 



[1] P. Braun-Munzinger and J. Stachel, Nucl. Phys. A638, 
3c (1998). 

[2] P. Danielewicz, Nucl. Phys. A661, 82c (1999). 
[3] E. Shuryak and O.V. Zhirov, Phys. Lett. B 89, 253 
(1979). 

[4] L. van Hove, Z. Phys. C 21, 93 (1983). 

[5] CM. Hung and E.V. Shuryak, Phys. Rev. Lett. 75, 4003 

(1995); Phys. Rev C 57, 1891 (1998). 
[6] D.H. Rischke, Y. Pursiin, J. A. Maruhn, H. Stocker, and 

W. Greiner, Heavy Ion Phys. 1, 309 (1996); D.H. Rischke, 

Nucl. Phys. A610, 88c (1996). 
[7] H. Liu for the E895 Collaboration, Nucl. Phys. A638, 

451c (1998). 

[8] I.N. Mishustin, V.N. Russkikh, and L.M. Satarov, Nucl. 



6 



Phys. A494, 595 (1989); Yad. Fiz. 54, 429 (1991) (trans- 
lated as Sov. J. Nucl. Phys. 54, 260 (1991)). 
[9] E.G. Nikonov, A. A. Sh anenko, and V.D Toneev, Heavy 
Ion Phys. 8, 89 (1998); [nucl-th/98020li . 
[10] V.D. Toneev, E.G. Nikonov, and A. A. Shanenko, in Nu- 
clear Matter in Different Phases and Transitions, eds. 
J. -P. Blaizot, X. Campi, and M. Ploszajczak, Kluwer 
Academic Publishers (1999), pp.309-320; Preprint GSI 
98-30, Darmstadt, 1998. 
[11] See e.g. F. Karsch, Talk given at the Strong and Elec 
troweak Matter '98, December 1998, Copenhagen, Eep- 



lat/9903031 



J. Zimanyi et al, Nucl. Phys. A484, 647 (1988). 
K. Redlich and H. Satz, Phys. Rev. D 33, 3747 (1986). 
C. Bernard et al, Nucl. Phys. (Proc. Suppl.) B47, 499 
(1996); ibid 503. 

L.M. Satarov, Sov. J. Nucl. Phys. 52, 264 (1990). 

J. Barrette et al, Phys. Rev. C 56, 3254 (1997). 

W. Reisdorf and H.G. Ritter, Ann. Rev. Nucl. Part. Sci. 

47, 663 (1997). 

N.N. Ajitanand et al, Nucl. Phys. A638, 451c (1998). 
N. Herrmann, FOPI Collaboration, Nucl. Phys. A610, 
49c (1996). 

J. Chance et al (The EOS Collaboration), Phys. Rev. 
Lett. 78, 2535 (1997). 

P.K. Sahu, W. Cassing, U. Mosel, and A. Ohnishi, Nucl. 
Phys. A672, 376 (2000). 

H. Appelshauser et al. (NA49 Collaboration), Phys. Rev. 
Lett. 80, 4136 (1998). 

S.A. Voloshin, Phys. Rev. C 55, 1632 (1997). 

S. S off, S.A. Bass, M. Bleicher, H. Stocker, and W. Grei- 



ner, micl-th/9903061 



J. Brachmann et al, Phys. Rev. C 61, 024909 (2000). 
J. Brachmann et al, Nucl. Phys. A619, 391 (1997). 
E.V. Shuryak, Phys. Rev. D 60, 115014 (1999); 
D. Teaney and E.V. Shuryak, Phys. Rev. Lett. 83, 4951 
(1999). 

P.F. Kolb, J. Sollfrank, and U. Heinz, Phys. Lett. B459, 
667 (1999). 

A.M. Poskanzer and S.A. Voloshin, Phys. Rev. C 58, 
1671 (1998). 



7 



